Low measurements of the proton form factor ratio ii^GeIGm 
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We present an updated extraction of the proton electromagnetic form factor ratio, /j^pGe/Gm , 
at low Q^. The form factors are sensitive to the spatial distribution of the proton, and precise 
measurements can be used to constrain models of the proton. An improved selection of the elastic 
events and reduced background contributions yielded a small systematic reduction in the ratio 
IJLpGe/Gm compared to the original analysis. 

PACS numbers: 13.0.Gp, 13.60.Fz, 13.88.+e, 14.20.Dh 
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I. INTRODUCTION /^pGe/Gm initially presented in Ref. pL], with improved 

We present a detailed reanalysis of polarization trans- 
fer measurements of the proton form factor ratio 
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selection of elastic events and significantly reduced con- 
tamination from quasielastic events in the target win- 
dows. The new results are typically lower by '^1%. 

The electric and magnetic form factors, Ge{Q'^) and 
Gm{Q'^)i describe the distribution of charge and magne- 
tization in the proton. The form factors are extracted in 
elastic electron-proton scattering and mapped out as a 
function of the four- momentum transfer squared, Q^, to 
yield the momentum-space structure of the proton. Pre- 
cision measurements of proton form factors over a large 
kinematic range can provide important constraints on 
models of the proton. However, when extracting the form 
factors from unpolarized cross section measurements us- 
ing the Rosenbluth separation technique, it is difficult to 
precisely separate Ge from Gm in the proton for very 
high or very low values. The addition of polariza- 
tion measurements allows for a much better sepa- 
ration oi Ge and Gm- Initial measurements for the pro- 
ton focused on the high-Q region which showed 
a significant falloff in the ratio /j^pGe/Gm with Q^, in 
contrast to previous extractions from Rosenbluth sepa- 
rations [11]. This difference is now believed to be due 
to the contribution of two-photon exchange effects which 
have a large impact on the extractions from the unpo- 
larized cross section measurements but have less impact 
on the polarization measurements [l2|-^| . These signifi- 
cantly improved measurements of G^; led to a great deal 
of theoretical work aimed at understanding this behav- 
ior [isl-llll, which showed, among other things, the im- 
portance of quark orbital angular momentum in under- 
standing the proton structure at high momentum [22142^ . 
These results also had a significant impact on studies of 
the correlations between the spatial distribution of the 
quarks and the spin or momentum they carry, showing 
that the spherically symmetric proton is formed from a 
rich collection of complex overlapping structures . 

While initial investigations focused on extending pro- 
ton measurements to higher Q^, the polarization mea- 
surements can also be used to improve extractions at low 
values, providing improved precision and less sensi- 
tivity to two-photon exchange corrections. The low-Q^ 
form factors relate to large-scale structures in the pro- 
ton's charge and magnetization distributions. As such, 
it has long been believed that the "pion cloud" contri- 
butions, e.g. the fluctuation of a proton into a virtual 
neutron-7r+ system, will be important at low Q^, as the 
mass difference means that the pion will contribute to 
the large distance distribution in the bound nucleon-pion 
system. It was recently suggested that such structures 
are present in all the nucleon form factors [2£], centered 
at ^ 0.3 GeV^, and that these structures reflect con- 
tributions from the pion cloud of the nucleon. However, 
the significance of the proposed structures and their in- 
terpretation as a pion cloud effect have been much dis- 
puted. This low-Q^ region is also important in parity vi- 
olating electron scattering measurements j27l - l3Q| aimed 
at investigating the strange-quark contributions to the 
proton electromagnetic structure. Isolating the strange 



quark contributions relies on precise determinations of 
the proton form factors at low Q^, including the impact 
of two-photon exchange corrections [31] (discussed fur- 
ther below). 

While the form factors encode information on the spa- 
tial structure of the proton, there are theoretical issues 
in extracting the spatial charge and magnetization distri- 
butions, discussed in detail elsewhere [32-36]. However, 
the difficulty in extracting true rest-frame distributions 
for the proton does not interfere with the comparison 
of form factor measurements and proton size/structure 
corrections to atomic levels in hydrogen. Extractions of 
the proton charge radius (37l-[4ll| define the proton root- 
mean-square (RMS) radius as the slope of the form fac- 
tor at = 0. This definition is consistent with the 
RMS radius needed in Lamb shift measurements in hy- 
drogen [42] and muonic hydrogen [43[. Corrections to 
the hyperfine splitting (ij-liG*! are also extracted directly 
from the form factors. The charge radius is of particular 
interest at present, due to the conflicting results between 
Lamb shift measurements on muonic hydrogen ^3[ and 
the electron scattering results and measurements from 
the Lamb shift in electronic hydrogen [42|. 

This experiment was motivated by the ideas discussed 
above: mapping out the large scale proton structure, the 
benefit of improved precision in proton form factors in or- 
der to extract strange quark form factors (and ultimately 
the proton weak radius) from parity violating measure- 
ments, and the importance of reducing the uncertainty in 
hyperfine splitting calculations arising from proton finite- 
size corrections. 



II. PREVIOUS MEASUREMENTS 

Since the 1960s, measurements of the unpolarized 
cross section for elastic e-p scattering have been used 
to separate Ge and Gm- The cross section is pro- 
portional to {rGj^ + ^G\)^ where r = Q'^ /Am^^ and 
£ = (1 + 2(1 + (QV4m2)) tan^ 0/2)-^. Keeping fixed 
while varying e allows for a "Rosenbluth separation" jl^ 
of the contributions from Ge and Gm- At high Q^, 
the factor of tG\^ dominates, as r becomes large and 
Gm > G\ (with Gm/Ge = /J^p at = q). This 
makes extraction of G^; difficult, as it contributes only 
a small, angle-dependent correction to the larger cross 
section contribution from Gm- Similarly, in the limit of 
very small Q^, and thus very small r, it is difficult to iso- 
late Gm except in the limit where £ ^ 0, i.e. scattering 
angle 180°. 

Polarization measurements are sensitive to the ratio 
Ge/Gm and thus, when combined with cross section 
measurements, can cleanly separate the electric and mag- 
netic form factors, no matter how small their contribu- 
tion to the cross section becomes. It has been known for 
some time that measurements of polarization ob- 

servables would provide a powerful alternative to Rosen- 
bluth separation measurements, but only in the last 
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decade or so have the high polarization, high intensity 
electron beams been available, combined with polarized 
nucleon targets or high efficiency nucleon recoil polarime- 
ters 

The ffist such measurements for the proton [gl, 
showed a decrease in i^pGe/Gm with Q^, which dif- 
fered from the existing Rosenbluth separation measure- 
ments, which showed approximate form factor scaling, 
i.e. I^pGe/Gm ^ 1. This discrepancy appeared to be 
larger than could be explained even accounting for the 
scatter in the previous Rosenbluth measurements 
A measurement using a modified Rosenbluth extraction 
technique [isj was able to extract the ratio /j^pGe/Gm 
with precision comparable to the polarization measure- 
ments, and showed a clear discrepancy, well outside 
of the experimental systematics for either technique. 
Experiments extending polarization measurements to 
higher show a continued decrease of i^pGe/Gm with 

It was suggested that the two-photon exchange (TPE) 
correction may be able to explain the discrepancy be- 
tween the two techniques [l2|, [l6|. While these correc- 
tions are expected to be of order aEM ^1%, they can 
have a similar e dependence to the contribution from Ge- 
Because the contribution to Ge is small at large Q^, a 
TPE correction of a few percent could still be significant 
in the extraction of Ge- It was estimated that a TPE 
contribution of ^5%, with a linear e dependence, could 
explain the difference [l2|, |49| , and early calculations sug- 
gested effects of a few percent, with just such a lin- 
ear ^-dependence [H, . These corrections should also 
modify the polarized cross section measurements, but it 
should be a percent-level correction in the extraction of 
Ge/ Gm , as there is is no equivalent amplification of the 
effect. Including the best hadronic calculations available 
yields consistency between the two techniques, and good 
separation oiGE and Gm up to high lH, [HI- Com- 
parisons of electron-proton and positron-proton scatter- 
ing can be used to isolate TPE contributions [51], and 
a series of such measurements are currently planned or 
underway [52l-[53|. 

At low values, the TPE should be well described 
by the hadronic calculations jisl, [55| , and in fact the con- 
tributions are smah for 0.3 < < 0.7 GeV^. While 
this is a region where high precision Rosenbluth separa- 
tions are possible, measurements prior to 2010 had rel- 
atively large uncertainties, typically 3-5% or more on 
I^pGe/Gm' Measurements using polarization observ- 
ables in this region can provide a significant improvement 
in precision, even in this low regime. The MIT-Bates 
BLAST experiment made measurements of /j^pGe/Gm 
usmg a polarized target [56^ for 0.15 < < 0.6 GeV^ 
with typical uncertainties around 2%, about a factor of 
two improvement over most earlier data. The experi- 
ment, which provided the best knowledge of the low 
proton form factor ratio when published, measured values 
below unity for > 0.2 GeV^, but concluded that the 
overall results were consistent with unity over the range 



of the experiment. Combined with the high-Q^ JLab 
data, which showed a clear deviation of i^pGe/Gm for 
^ 0.8 GeV^ , this suggested that the ratio was unity 
at very low and then began to fall somewhere in the 
range of 0.2-0.7 GeV^. The fact that there was no clear 
indication of where the ratio began to fall below unity 
was one of the motivating factors for this measurement. 
The updated results of this reanalysis of [1] provide an 
independent extraction of /j^pGe/Gm in this kinematic 
region, with precision comparable to the BLAST results. 
More recently, JLab experiment E08-007 [4l|, a high- 
statistics follow up to the work we present here, used 
the same polarization transfer techniques but with coin- 
cident detection of the final-state electron and proton for 
all kinematics, yielding measurements of i^pGe/Gm with 
average uncertainties below 1.2%. 

Last year, new measurements in this region were 
also obtained by an experiment at Mainz [40]. The 
experiment made high-precision measurements of un- 
polarized cross sections at ~1400 kinematic points for 
< 1 GeV^. While they do not provide direct Rosen- 
bluth extractions of G^; and Gm, they show a global fit 
to their cross section results. Their extraction of Gm is 
systematically 2-4% above previous world's data, imply- 
ing a difference of 4-8% in the extrapolation of the cross 
section to e = 0. It is difficult to determine how much 
of their error band could be strongly correlated in Q^, 
as there is no information given on the size or sources of 
systematic uncertainty assumed in their analysis. While 
they apply a very limited form of the two-photon ex- 
change corrections [57], which is neglected in most pre- 
vious extractions, this should only reduce their value of 
Gm relative to the uncorrected results, implying that the 
true discrepancy is even larger. At this point, it is not 
clear why there is such a large discrepancy between their 
fit and previous measurements. 



III. EXPERIMENT DETAILS 

This experiment was carried out in Hall A of the 
Thomas Jefferson National Accelerator Facility (JLab), 
in the summer and fall of 2006, as part of experiment 
E05-103 [58j. While the experiment was focused on po- 
larization observables in low energy deuteron photodis- 
integration [59], elastic electron-proton scattering mea- 
surements used to calibrate the focal plane polarimeter 
provided high statistics data that allowed for an improved 
extraction of the proton form factor ratio /j^pGe/Gm at 
low 

A polarized electron beam was incident on a cryogenic 
liquid hydrogen target, nominally 10 cm in length for 
the 362 MeV beam energy running and 15 cm for the 
687 MeV settings (the target length was misstated as 
15 cm for all runs in the previous publication [1]). The 
target cells are Al, with beam entrance windows about 
0.1 mm thick, and beam exit and sides ~0.2mm thick 
(with some variation between the different targets). Elas- 
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tic e-p scattering events were identified by detecting the 
struck proton in one of the High Resolution Spectrom- 
eters (HRS) [60]. Data were taken with a longitudinal 
polarization of approximately 40% and with the beam 
helicity flipped pseudo-randomly at 30Hz. For some set- 
tings, the scattered electron was detected in the other 
HRS spectrometer. 

The polarization of the struck protons is measured in 
a focal plane polarimeter (FPP) in the proton spectrom- 
eter. Operation and analysis of events in the FPP is 
described in detail in Refs. [1, [60|. Analysis of the angu- 
lar distribution of rescattering in the polarimeter allows 
us to extract the transverse polarization at the detector, 
which can be used to reconstruct the longitudinal and 
transverse (in-plane) components of the polarization of 
the elastically scattered protons. In the Born approx- 
imation, the ratio of these polarization components is 
directly related to the ratio Ge/Gm^ 

where Cz^x are the longitudinal and transverse compo- 
nents of the proton polarization, Eq is the beam energy, 
and Oe and E' are the scattered electron's angle and mo- 
mentum (reconstructed from the measured proton kine- 
matics). Because the extraction of i^pGe/Gm depends 
on the ratio of two polarization components, knowledge 
of the absolute beam polarization and FPP analyzing 
power are not necessary, although high polarization and 
analyzing power improve the figure of merit of the mea- 
surement. 

In the experiment, we measure the polarization not at 
the target, but in the spectrometer focal plane, and the 
asymmetry in the rescattering is sensitive only to polar- 
ization components perpendicular to the proton direc- 
tion. If we look at the central proton trajectory, where 
the spectrometer is well represented by a simple dipole, 
then the transverse component, Cx^ will be unchanged, 
while the longitudinal component, Cz^ will be precessed 
in the dipole field. If we chose a spin precession angle, 
X, near 90 degrees, the longitudinal and transverse polar- 
ization components at the target will yield "vertical" and 
"horizontal" components in the frame of the focal plane 
polarimeter, allowing for both to be extracted by a mea- 
surement of the azimuthal distribution of rescattering in 
the carbon analyzer. In the analysis, we use a detailed 
model of the spectrometer to perform the full spin pre- 
cession, rather than taking a dipole approximation, as 
described in detail in Ref. [8]. 

A follow-up experiment, JLab E08-007 [61] was pro- 
posed to make extremely high precision measurements in 
this kinematic regime. The measurement was run in the 
summer of 2008, and in the analysis of the E08-007 data, 
it was observed that the result was somewhat sensitive to 
the cuts applied to the proton kinematics when isolating 
elastic e-p scattering. 

In this experiment, only the proton was detected for 
most kinematic settings, and the elastic scattering events 



TABLE I: Kinematics and FPP parameters for the measure- 
ments. Oj^^ and Tp are the proton lab angle and proton ki- 
netic energy, respectively. Tanaiyzer is the thickness of the 
FPP carbon analyzer and x is the spin precession angle for 
the central trajectory. The final column shows which kine- 
matics had single-arm (S), coincidence (C), or a combination 
of both (C/S). 
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were isolated using cuts on the over-determined elastic 
kinematics. In the original analysis [1], relatively loose 
cuts were applied because the measurement was statis- 
tics limited and little cut dependence had been observed 
in previous measurements 0, [sl, |62| . Most of these mea- 
surements had high resolution reconstruction of both the 
proton and electron kinematics, and so loose cuts on the 
combined proton and electron kinematics provided clean 
isolation of the elastic peak. In addition, the previous 
measurements were generally at higher and so of sig- 
nificantly lower statistical precision, typically 3-5%, so it 
was difficult to make precise evaluation of the impact of 
tight cuts on the proton kinematics. Because the elastic 
events could be cleanly identified without tight cuts on 
the proton kinematics, this was not considered to be a 
significant concern. 

In the follow-up experiment, E08-007 [61], the electron 
was detected in a large acceptance spectrometer with lim- 
ited momentum and angle resolution. The electron detec- 
tion led to significant suppression of scattering from the 
target windows, but the poor electron resolution required 
that the elastic peak be defined using cuts on the proton 
kinematics. Because of this, and the high statistics of 
the data set, it was possible to make detailed studies of 
the cut dependence of the result. It was found that there 
were small but noticeable changes in the extracted form 
factor ratio if the proton kinematic cuts were made too 
loose, even in cases where the endcap contributions were 
small. 

Motivated by these issues, we reanalyzed the data from 
our experiment. We include a more careful examination 
of cuts used to identify the proton events and an updated 
evaluation of the contribution from the target endcaps. 
With our new, more restrictive cuts, there were small but 
systematic changes in the extracted form factors. These 
were mainly due to the reduction in the contribution from 
electron scattering in the Aluminum endcaps rather than 
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any changes in the events corresponding to scattering 
from hydrogen. 

One of the most important issues in the original analy- 
sis was the correction for events that came from scatter- 
ing in the Aluminum endcaps of the targets, and there 
were several difficulties involved in making these correc- 
tions. For systematic checks, we took data with the elas- 
tic peak centered on the focal plane, but also with the 
spectrometer momentum approximately 2% higher and 
lower, to map out the response of the FPP across the 
focal plane. Data on the Aluminum dummy targets were 
typically taken for only one setting, and so there was 
a systematic uncertainty associated with the stability of 
the size of the background contamination and the possi- 
ble variation of the polarization transfer coefficients mea- 
sured from the dummy target. Data were taken using 
both 15 cm and 10 cm cryogenic hydrogen targets, but 
only dummy foils for the standard 4 cm and 15 cm cry- 
otargets were available, yielding some additional system- 
atic uncertainties for the 10 cm targets. Finally, for some 
runs the beam position was not perfectly stabilized on the 
dummy foils and the beam, rastered to a 4 by 4 mm^ spot 
at the target, either partially missed the dummy foils or 
impinged on both the 4 cm and 15 cm dummy targets. 
While the beam position is continuously monitored and 
we correct for any deviation in the event reconstruction 
based on the position, the luminosity is not well known 
if the beam is partially missing the foils. Therefore, the 
relative normalization of the contribution from the target 
endcaps and the dummy foils had to be determined by 
looking at quasielastic events that are above the thresh- 
old for scattering from the proton, rather than being cal- 
culated directly, yielding an additional systematic in the 
relative normalization of the endcaps and dummy foils. 

In the original analysis, endcap scattering typically 
yielded 3-5% of the cross section after all cuts were ap- 
plied (much less for the two coincidence settings), so there 
was a small but significant correction. Because of the is- 
sues mentioned above, there were very large systematic 
uncertainties associated with these corrections. In ad- 
dition, the original analysis applied the full set of cuts 
for elastic scattering to data from the dummy target, 
yielding measurements of the polarization coefficient for 
endcap scattering with extremely poor statistics and thus 
large fluctuations. We now use much tighter cuts on the 
reconstructed target position to try to remove most of 
the endcap contributions, resulting in contributions of 
^0.5%. While the cuts reduce the statistics of the main 
measurement somewhat, the final uncertainty is often 
better, as the endcap subtraction, which had large statis- 
tical and systematic uncertainties, is now much smaller. 
We also use looser cuts when extracting Cx and Cz from 
scattering in the aluminum endcaps, with an extra sys- 
tematic uncertainty applied to account for possible cut 
dependence. 



IV. ANALYSIS DETAILS 

For elastic scattering using an electron beam with a 
known energy the complete scattering kinematics can be 
determined from the measurement of a single kinematic 
quantity, typically the angle or energy of the final state 
electron or proton. If two quantities are measured, then 
the consistency of the two kinematic variables can be 
used to determine if the event was associated with elastic 
scattering. For this analysis, we use the proton scattering 
angle and momentum to reconstruct the kinematics and 
to identify elastic events. For some kinematics, the elec- 
tron was also detected, which allows for almost complete 
suppression of events coming from quasielastic scattering 
in the aluminum entrance and exit windows of the target. 
To identify the elastic peak, we use the difference between 
the measured proton momentum and the momentum cal- 
culated based on the measured proton angle measured 
proton angle. The specific variable we use is DpKin, 
which is the momentum difference, Pp — Peiasud^p), di- 
vided by the central momentum setting of the proton 
spectrometer. This yields a fractional momentum devia- 
tion from the expectation for elastic scattering. 




-0.06 -0.04 -0.02 0. 0.02 0.04 0.06 



FIG. 1: Reconstructed target position Ytg vs. DpKin = 
{Pp — Peiastic{Op)) /phrs for the measurements on the 10 cm 
liquid hydrogen target (top) and on the 4 cm and 15 cm Alu- 
minum "dummy" foils (bottom). Note that Ytg is the position 
transverse to the spectrometer optic axis, not the position 
along the beamline; this difference leads to the target dimen- 
sions being reduced by a factor of ^2 here. The elastic peak 
is clearly visible at DpKin ~ for the LH2 target, while the 
broad quasielastic contributions from endcap scattering are 
visible at the ends of the LH2 target. 
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Figure [T] shows the distribution of events versus DpKin 
and the reconstructed target position, 1^^, as seen by the 
spectrometer. For the hydrogen target (top panel), there 
is a strong peak at DpKin ~ 0, corresponding to elastic 
events. At the extreme Ytg values, there is a faint but 
broad distribution corresponding to quasielastic scatter- 
ing in the endcaps. We apply a cut to Ytg to remove most 
of the contribution from the endcap scattering, and use 
the measurements from the 4 cm and 15 cm dummy tar- 
get (bottom panel) to subtract the residual contribution. 
Note that for the spectra shown in Fig. (TJ the length of 
the LH2 target does not match either the inner or the 
outer pair of foils from the dummy target. This means 
that the acceptance as a function of DpKin depends on 
Ytg and so will not be identical for the endcaps and the 
foils in the dummy target. This is clear for the outer foils 
of the dummy target, where there is a significant loss of 
events at extreme positive (negative) values of DpKin 
for the upstream (downstream) dummy foils. 

For each setting, three measurements were taken; 
one with the elastic peak positioned at the central mo- 
mentum of the spectrometer, and two where the elastic 
peak was shifted up (down) by 2% in momentum. This 
allowed us to verify that the result was independent of 
the position of the events on the focal plane. However, 
dummy events were typically taken at only one of these 
three settings, and the extracted endcap contribution and 
quasielastic recoil polarizations taken from that measure- 
ment were applied to all three settings, so the DpKin 
distributions will not be exactly identical, especially far 
away from the elastic peak. The dummy spectra are nor- 
malized to match the observed "super-elastic" contribu- 
tion {DpKin > 0.03 in Fig. [2]) in the LH2 data, using 
only the inner foils for the dummy target, as they have 
a DpKin acceptance which better matches the endcaps. 
Figure [2] shows the spectra for the LH2 target (thick black 
histogram) and the dummy target (grey histogram), af- 
ter the dummy target has been normalized in the region 
indicated by the vertical dashed lines. After normalizing 
the spectra in this region, we can determine the endcap 
contribution under the elastic peak. The region used to 
define elastic events in the analysis is indicated by the 
vertical dotted lines. We take a conservative approach 
and apply a 50% systematic uncertainty to the size of 
the endcap contribution when making the correction for 
these events to account for the impact of the different 
DpKin spectra between the endcaps and the dummy foils 
and possible variation for the settings which are shifted 
by ±2% in momentum. 

Having determined the contribution from endcap scat- 
tering, we use the data from the dummy targets to de- 
termine the contributions from quasielastic scattering to 
the recoil polarization components Cx and Cz- If we ap- 
ply the same cuts to the dummy target as we use in the 
analysis of the hydrogen, there is very little data left, and 
we can not make a reliable extraction of Cx and Cz - For 
the quasielastic scattering, we use all four aluminum foils 
and a broader cut on DpKin to determine the quasielas- 
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FIG. 2: (Color online) The DpKin distribution for the hy- 
drogen target (thick black histogram) and the dummy targets 
(thin red histogram). The dashed vertical lines indicate the 
region used to normalize the dummy contribution to match 
the contribution from the aluminum endcaps of the hydrogen 
target and the vertical dotted lines indicate the part of the 
elastic peak used in the analysis. 



tic values for Cx and C^, and then assume that the co- 
efficients are identical when looking at the central part 
of the quasielastic spectrum. Comparisons showed com- 
plete consistency between the extracted values of Cx and 
Cz when comparing the inner and outer dummy foils for 
all kinematics, or when comparing the central part of 
the quasielastic peak to the off-peak contributions. Be- 
cause we could not make a precise determination of Cx 
and Cz without averaging over a larger kinematic region, 
we apply an uncertainty to Cx and Cz of 0.02 and 0.05 
respectively, compared to typical values for these polar- 
ization components in this experiment of 0.08-0.2 for Cx 
and 0.15-0.3 for C^. 

In the original analysis [1], the Ytg cut was loose and so 
there was a large (3-5%) contribution from endcap scat- 
tering which had to be subtracted. Because the tight cuts 
used on the elastic events were also applied to the dummy 
spectra used to subtract endcap scattering contributions, 
the statistical uncertainty on these subtractions could be 
very large. Therefore, ffuctuations in the low statistics 
dummy measurements led to large uncertainties and sig- 
nificant fiuctuations in the dummy-subtracted measure- 
ments. In the present analysis, the endcap contributions 
are greatly reduced, with a maximum contribution well 
below 1%, such that the conservative systematic uncer- 
tainties assumed for the dummy normalization and po- 
larization coefficients yield only small uncertainties in the 
final result. While the tighter cuts yield slightly reduced 
statistics in the elastic peak, the total statistical uncer- 
tainty is sometimes smaller because the background con- 
tribution was reduced. Note that for a few settings, ad- 
ditional runs were included, improving the statistics by 
5-15%, but this was a small effect compared to the mod- 
ified cuts. 

There were also some small changes in the evaluation 
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of the systematic uncertainties. In the previous analysis, 
the systematic uncertainty from the endcap contribution 
was folded into the reported statistical uncertainties, and 
these are now part of the quoted systematics. In addi- 
tion, the estimated systematics are somewhat larger than 
in the previous analysis, due to a more detailed analy- 
sis of the uncertainty in the spin precession through the 
spectrometer [63]. 

The proton energy loss, which can be significant for 
the low kinematics, was also more carefully evalu- 
ated, leading to a small change in the average for 
each bin. For the 362 MeV running, where the proton 
was detected at small angles, the energy loss depends on 
the position in the target where the scattering occurs. 
Figure [3] shows DpKin vs. Ytg for one of the 362 MeV 
runs with an average proton energy loss is applied to all 
events. Positive Ytg values correspond to the upstream 
portion of the target, where all events exit through the 
side of the target and travel through a constant amount 
of hydrogen and aluminum and can be well corrected 
assuming a fixed energy loss. Events that exit through 
the downstream end of the target lose less energy because 
they pass through less material, yielding a F^^-dependent 
position for the elastic peak. This yields a reduced proton 
energy loss, and thus a higher apparent proton momen- 
tum, for events that occur near the exit window. For 
the kinematics where this effect is important, we apply 
a F^^-dependent cut, cut corresponding to a two-sigma 
region around the elastic peak for each region of Ytg^ as 
indicated by the graphical cut displayed in Fig.[3l The re- 
constructed value of DpKin is only used to select elastic 
events, so while a position-dependent energy loss could 
have been applied, one would still end up with the same 
set of good events passing the cuts. In our approach, we 
are not sensitive to any imperfections in the energy loss 
correction, since we use a two-sigma cut for all Ytg values. 



TABLE II: Kinematic-dependent cuts applied to the data. 
The Ytg cut is chosen to significantly suppress any contribu- 
tions from the target endcap (as shown in Fig. [3|). 





^ 

^lab 


Ytg cut 


6p/p cut 


(GeV^) 


(deg) 


(cm) 




0.215 


28.3 


-0.022 < Ytg < 0.018 


\6p/p\ < 0.045 


0.235 


23.9 


-0.022 < Ytg < 0.018 


\6p/p\ < 0.045 


0.251 


18.8 


-0.018 < Ytg < 0.012 


\6p/p\ < 0.045 


0.265 


14.1 


-0.014 < Ytg < 0.010 


\6p/p\ < 0.045 


0.308 


47.0 


-0.025 < Ytg < 0.020 


\6p/p\ < 0.040 


0.346 


44.2 


-0.025 < Ytg < 0.020 


\6p/p\ < 0.040 


0.400 


40.0 


-0.028 < Ytg < 0.022 


\6p/p\ < 0.040 


0.474 


34.4 


-0.024 < Ytg < 0.020 


\6p/p\ < 0.040 




6000 



The kinematic-dependent cuts are detailed in Table ITTl 
In addition, several cuts were applied to all kinemat- 
ics. A cut was applied on the out-of-plane angle, \Otg\ < 
0.06 rad, and the in-plane angle, < 0.03 rad, to 

ensure events were inside of the angular acceptance of 



).06 -0.04 -0.02 0. 0.02 0.04 0.06 
DpKin 

FIG. 3: (Color online) The reconstructed target position, Ytg 
vs. DpKin, the deviation of the momentum from the ex- 
pected elastic peak position. A correction for the average 
energy loss is applied, but there is a significant difference for 
events on the upstream side of the target, which exit through 
the side wall of the target, and events which occur nearer 
the downstream end of the target and have less energy loss. 
The band indicates the graphical cut placed on these runs, to 
approximate a two-sigma range for each Ytg value. 



the spectrometer. A two-sigma cut was applied on the 
DpKin peak, with a Ft^-dependence cut for the low en- 
ergy kinematics to account for the position-dependent 
average energy loss as shown in Figure [3l The tracks 
before and after the Carbon analyzer were used to deter- 
mine the scattering location and scattering angle in the 
analyzer. Events were required to have the secondary 
scattering occur within the analyzer, and angle between 
5 and 50° were accepted. In addition, we apply a cone 
test [8] to ensure that there is complete azimuthal ac- 
ceptance in the EPF. We do this by requiring that the 
EPF would have accepted events with any azimuthal an- 
gle given the reconstructed vertex and scattering angle. 
This ensures that any asymmetry in the acceptance or 
distribution of events does not lead to a difference in the 
scattering angle distribution for vertical and horizontal 
rescattering. A significant difference between the rescat- 
tering distribution for vertical and horizontal rescattering 
events would yield a different average analyzing power, 
and the analyzing power would not cancel out in the ratio 
of polarization components. 

The combination of the more restrictive cuts on the 
elastic events and the associated reduction in contami- 
nation due to scattering from the target windows leads 
to a reduction in the extracted ratio that is typically 
at the 1-2% level. The largest effect is due to the im- 
proved correction for endcap scattering, mainly due to 
cuts which significantly reduced the size of this contri- 
bution. There is also a 1% reduction in the coincidence 
settings, where there are negligible endcap contributions, 
which is due to the tighter cuts on the proton kinemat- 
ics. Tight elastic kinematics cuts using just the proton 
will remove events where there is a larger than average 
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TABLE III: Systematic uncertainties on i? = /j^pGe/Gm ■ See 
text for details. 





SR 


5R 


SR 


V ; 


( pnrlr^^D 1 


( DDtirs ) 


( c^^t^\ 


0.215 


0.0012 


0.0079 


0.0141 


0.235 


0.0004 


0.0079 


0.0120 


0.251 


0.0003 


0.0078 


0.0107 


0.265 


0.0003 


0.0076 


0.0098 


0.308 




0.0091 


0.0077 


0.346 




0.0086 


0.0066 


0.400 


0.0010 


0.0088 


0.0056 


0.474 


0.0007 


0.0117 


0.0049 



error in the reconstruction of the proton scattering an- 
gle or momentum due to multiple scattering or imperfect 
track reconstruction. While these errors are small, the 
reconstructed kinematics are used to determine the spin 
propagation through the spectrometer, and thus the im- 
pact of the poor reconstruction may be amplified in eval- 
uating the spin precession. Table Hill shows the various 
contributions to the systematic uncertainty as a function 
of Q^. At high Q^, the uncertainty in the spin precession 
due to imperfect knowledge of the spectrometer optics 
dominates. At low Q^, the uncertainty is dominated by 
our ability to determine the cut dependence of the re- 
sult. The cut-dependent uncertainties come mainly from 
two sources; possible variation of the result due to the 
cuts on ytg and DpKin. While no systematic cut depen- 
dence with the ytg cut was observed, we apply a 0.4% 
uncertainty as a conservative estimate based on examin- 
ing the variation of /j^pGe/Gm with the ytg cut, in par- 
ticular for the coincidence data where the background 
contributions are smaller. For DpKin, we estimate the 
uncertainty based on varying the width of the cut around 
the elastic peak. This was done for both these data and 
the E08-007 results and no systematic cut dependence 
was observed, so the scatter of the results was taken as 
a conservative estimate of the systematic uncertainties. 
Because the data taken at low beam energy do not have 
sufficient statistics to set precise limits, we fit the un- 
certainties from the higher energy measurements and the 
E08-007 results and find a behavior consistent with l/Q^, 
which we use to obtain the quoted uncertainties for the 
low values. 



V. RESULTS 

The results of the reanalysis are given in Table IIVI and 
shown in Figure HI which presents the updated results 
along with previous measurements and a selection of fits. 
The updated analysis yielded a systematic decrease of 
^1% in the extracted ratio, except for the highest 
point which decreased by 5%. The analyzing power has 
been extracted from these data [64], but the quality of 



this extraction does not impact these results, as the ana- 
lyzing power cancels out in the ratio of Eq. [TJ However, 
because the FPP efficiency and analyzing power are sig- 
nificantly lower for the data taken at a beam energy of 
362 MeV (due to the lower proton momentum and thin- 
ner analyzer), the statistical uncertainty in these points is 
much larger. The two-photon exchange corrections from 
the hadronic calculation of Blunden, et al. [55] are 0.35% 
for the data below = o.3 GeV^ and 0.2% for the 
higher points. This is well below the statistical and 
systematic uncertainties for all points, and no correction 
(or uncertainty) for the TPE effects is included in the 
extraction. 



1.05 




Earlier measurements 
Punjabi et al. 
Crawford et al. 
Paolone et al. 
X. Zhan et al. 
This work 



FIG. 4: (Color online) The proton form factor ratio as a func- 
tion of (with the four low-Q^ measurements combined into 
one data point) shown with previous extractions with total 
uncertainties <3%. The curves are various fits [7j,ll3|,|26|,|4l|, 
while the dot-dashed curve and associated error bands show 
the result of the fit to the recent Mainz measurements [401]. 



The results in Fig. |4] show that the original conclu- 
sions of [1] are largely unaffected. The new results sup- 
port even more clearly the conclusion that the ratio 
IJ^pGe/Gm is below one even for these low values, 
with the change from previous Rosenbluth separations 
being driven mainly by a change in Ge-, with a smaller 
change in Gm- The previous hint of a local minimum 
near = 0.35 — 0.4 GeV^ was a consequence of the 
point near 0.5 GeV^, and there is no longer any indica- 
tion for this in our measurement. These results further 
support the observation that the decrease of the ratio 
below unity occurs at low Q^, and thus we expect that 
there will be a slightly larger impact on the extraction of 
strange quark contributions, as discussed in the original 
paper [1]. 

A comparison of the high-precision measurements at 
low-Q^ shows some small but systematic differences. The 
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TABLE IV: Experimental Results. R is given along with its 
statistical and systematic uncertainties. The last column (f) is 
the fractional contribution from scattering in the target end- 
caps, along with the statistical uncertainty; a 50% systematic 
uncertainty is also applied. The contribution is negligible for 
the coincidence settings. For = 0.474 GeV^, dummy mea- 
surements were taken at all three sub-settings, and the range 
of results is given. ^The final entry is the average of the four 
low- statistics point below = 0.3 GeV^. 





R — jj^pGe/Gm 


f 


(GeV") 




(%) 


0.215 


0.8250 lb 0.0483 ± 0.0162 


0.26(3) 


0.235 


0.9433 ± 0.0414 ± 0.0144 


0.13(2) 


0.251 


0.9882 zb 0.0420 ± 0.0132 


0.19(3) 


0.265 


0.9833 =b 0.0349 ± 0.0124 


0.16(2) 


0.308 


0.9320 =b 0.0123 ± 0.0119 




0.346 


0.9318 ± 0.0098 ± 0.0108 


-/0.40(2) 


0.400 


0.9172 ± 0.0109 ± 0.0105 


0.65(4) 


0.474 


0.9225 =b 0.0160 ± 0.0127 


0.4-0.6 


0.246^ 


0.9465 ± 0.0204 ± 0.0137 


n/a 



results from the Mainz cross section measurements [4Q| 
are 1-2% above the recoil polarization measurements 
from this work and the lower results from the re- 
cent JLab E08-007 measurement [41], although they are 
in agreement with the E08-007 results at higher val- 
ues. One concern for the results extracted from the 
Mainz cross section measurements is the sensitivity to 
two-photon exchange (TPE) corrections [15]. For the 
kinematics of the Mainz experiment, these corrections 
are fairly small, ^2%, but this is very large compared to 
the statistical (<0.2%) and systematic (<0.5%) uncer- 
tainties applied in the global fit to Ge and Gm- Thus, 
if ignored, this could yield significant corrections com- 
pared to the quoted uncertainties. Coulomb corrections 
were applied using the prescription of McKinley and Fes- 
hbach [53, which corresponds to the = limit of 
the Coulomb distortion correction (the soft-photon ap- 
proximation of the full TPE corrections). However, over 
much of the range of the experiment, applying the 

= correction is worse than neglecting the correction 
altogether, as the Coulomb correction changes sign at 

~ 0.15 GeV^ [65]. An estimate of the impact of these 
corrections [66] suggests that an improved prescription 
would decrease jj^pGe/Gm extracted from the cross sec- 
tions by 1-3% for ^ 0.1 GeV^ in a direct Rosenbluth 
separation, although a more complete analysis is required 
to determine the impact on their global fit. Bernauer et 
al.^ have examined the impact of these corrections [67] 
in more detail, although the TPE prescription they ap- 
ply @ is only valid up to = q.I GeV^. They find 
that their extracted value for /j^pGe/Gm changes by more 
than the total quoted uncertainty for all values up to 
0.1 GeV^, and while they find smaller changes at larger 

values, this is where the prescription is not expected 



to be valid. 

If the Coulomb corrections bring the Rosenbluth [40| 
extractions into agreement with the recoil polarization 
data, there is still a small systematic disagreement be- 
tween these and the polarized target measurements from 
BLAST [5^. At this point, we are unaware of any the- 
oretical argument that would explain a difference be- 
tween the results of the two different polarization tech- 
niques. This discrepancy can be further examined in 
the second phase of the JLab E08-007 experiment [6l[ 
which will make extremely high precision measurements 
of iIpGe/Gm down to ^ 0.015 GeV^, allowing for 
a comparison with the BLAST measurements using the 
same basic technique. 
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FIG. 5: (Color online) The proton form factor ratio 
/J^pGe /Gm vs. Q^, compared to several low-Q^ models. The 
curves shown are Miller's light-front cloudy-bag model calcu- 
lation 69]; Boffi's point-form chiral constituent quark model 
calculation 70]; Faessler's light-front quark model calcula- 
tion 71]; Lomon's vector- meson dominance model |7^; the 
dispersion analysis of Belushkin, Hammer and Meissner [73]; 
and the model of de Melo et al [tI. 



Figure [5] shows the measurements compared to a set of 
theoretical curves. The first type of calculation is based 
on the constituent quark models, which was quite suc- 
cessful in describing the ground state baryon static prop- 
erties. To calculate the form factors, relativistic effects 
need to be considered. Miller |69| performed a calcula- 
tion in the light-front dynamics including the effect from 
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the pion cloud. Boffi et al [70| performed a point form 
calculation in the Goldstone boson exchange model with 
point-like constituent quarks. Faessler et al [7l| used a 
chiral quark model where pions are included perturba- 
tively and dress the bare constituent quarks by mesons 
in a Lorentz covariant fashion. Another group of calcu- 
lations is based on the Vector Meson Dominance (VMD) 
picture, in which the scattering amplitude is written as 
an intrinsic form factor of a bare nucleon multiplied by 
an amplitude derived from the interaction between the 
virtual photon and a vector meson. This type of mod- 
els usually involve a number of free parameters for the 
meson mass and coupling strength. Lomon |75|, [76| 
performed the VMD fits by including additional vector 
mesons and pQCD constraints at large Q^. Belushkin et 
al. [73] performed a calculation using dispersion relation 
analysis with additional contribution from pir and KK 
continua. More recently, de Melo et al. [74] performed a 
calculation in the light-front VMD model by considering 
the no n- valence contribution of the nucleon state. While 
most of the theoretical curves are a few percent higher, 
the calculations of Miller [69| and de Melo et al. [74] gen- 
erally reproduce the large deviation from jUpGe/Gm = 1 
in this low region, emphasizing the pion cloud or non- 
valence effect. 

Our new results, as well as other high-precision mea- 
surements at low Q^, show that /j^pGe/Gm < 1 even 
down to very low values of Q^. A global fit [i^ to the 
cross section and polarization measurements in this re- 
gion, including the data presented in this work, indicates 
that G^; is ^2% below previous fits that did not include 
the low polarization measurements, while Gm is ap- 
proximately 1% higher. These small changes in the low 
form factors impact other measurements as well. For 
example, it was recently pointed out [77] that the re- 
duction in the form factor yields an agreement between 
studies of the asymmetry in the D(e,e'p)n reaction at low 
missing momentum in polarized target measurements at 
NIKHEF and MIT-Bates [78-80], and recent measure- 
ments at Jefferson Lab [81]. Similarly, this small shift in 
Ge and Gm at low modifies the expected asymme- 
try in parity- violating elastic electron-proton scattering, 
which serves as the baseline when extracting the strange- 
quark contribution to the proton form factors [27-30]. 
The effect is relatively small for any given extraction, 
especially at forward angles where there is a partial can- 
cellation due to the changes in G^; and Gm- However, 
because this is a systematic correction to all such mea- 
surements, the updated form factors could have a small 
net contribution on the extracted strange-quark contri- 
butions. 

The form factors at very low are important in both 
extracting the proton size and as input to other finite-size 
corrections in atom ic phys ics, e.g. the Zemach radius in 
hyperfine splitting [44-[46j. Measurements of the ratio 
IJ^pGe/Gm are not, by themselves, sufficient to extract 
the proton radius or calculate the finite-size corrections, 
as these require the individual form factors Ge{Q'^) and 



Gm{Q'^)' However, the polarization measurements can 
be used to improve global fits to cross section data, in 
which determining the relative normalization of differ- 
ent data sets is often the limiting factor in the system- 
atic uncertainties [37]. It is often possible to modify the 
normalization of certain data sets and modify the ratio 
IJ^pGe/Gm in such a way that the cross section mea- 
surements are still relatively well fit. Thus, having direct 
constraints on /j^pGe/Gm with high precision allows for a 
better determination of these normalization factors, and 
an improved extraction of the form factors. The fit of 
Ref. [13] was updated in Ref. [41] to include the data 
presented here as well as additional polarization mea- 
surements [U, [sil, yielding improved extractions of 
Ge^ Gm^ as well as the proton charge and magnetization 
radii. While these polarization measurements do not go 
as low in as one would like for the extraction of the 
radius, they nonetheless play an important role in the 
extraction, mainly by providing improved determinations 
of the relative normalization of the different experiments. 
In addition, higher order terms in the expansion need 
to be constrained to obtain a measure of the radius [s^] • 

This updated global analysis yielded an RMS charge 
radius of 0.875(10) fm, consistent with the CODATA 
value [42] and recent extraction from cross section mea- 
surements from Mainz [40]. Combined, these indepen- 
dent extractions based on the electron-proton interaction 
yield a radius of 0.8772(46) fm [4l|, more than seven stan- 
dard deviations from the recently pubhshed PSI muonic 
hydrogen Lamb shift measurement [43] of 0.8418(7) fm. 
While higher order corrections, which depend on the val- 
ues of the form factors at finite Q^, can modify the ex- 
tracted radius, these corrections appear to be far too 
small to explain the discrepancy, although this is still 
under investigation [ss'-^SG']. A recent work [ST^ has pro- 
posed a possible mechanism to explain a difference be- 
tween electronic and muonic probes of the proton struc- 
ture, due to off-shell effects in the hadronic intermediate 
state in the two-photon exchange diagrams. However, 
while this is an area that has received a great deal of at- 
tention in the recent past |83l-[85l, [88l - [92| , the question is 
still unresolved. 

The updated global fit of the low data yields a mag- 
netic RMS radius of 0.867(20) fm [4l| . This is in reason- 
able agreement with other extractions, 0.85(3) fm from 
Kelly [32] and 0.857 from Hammer and Meissner |93j. 
However, neither of these included Coulomb distortion 
corrections, and they are global fits that include a large 
body of high-Q^ data which do not contain useful in- 
formation on the proton radius, and thus may not be 
reliable when it comes to examining the radius. Sick ^3 
extracted the Zemach radius, which depends on both the 
charge and magnetic distributions, and converting this to 
a magnetic radius yields a value of 0.855(35) fm [93], in 
agreement with the updated fit. However, the data used 
by Sick are also included in this global fit, so these can- 
not be considered independent extractions. The Mainz 
result [io^ is 0.777(17) fm, weh below the other electron 
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scattering analyses, but the strong dependence of the 
Coulomb distortion at very low Q^, neglected in their ex- 
traction, may have a significant impact on the extraction 
of the radius [66]. An updated estimate yields a 1.5cr shift 
in their magnetic radius (and almost no change in the 
electric radius), yielding 0.803(17) fm [67], where the un- 
certainty does not include any contribution related to the 
two-photon exchange corrections. This increased radius 
is now only 2Aa from our result, but this updated value 
still appears to be an underestimate, as their weighted 
averaging of results from different fits [94] leads to a bias 
towards fits with fewer parameters. They give systemat- 
ically lower values for the magnetic radius than the fits 
with more parameters, as one might expect given that 
the precision of the extraction of Gm tends to improve 
as increases, and so inclusion of the high-Q^ data can 
dominate the fit if there is not enough flexibility in the 
fit. This does not appear to affect their extraction of the 
charge radius, as the extracted radius does not show the 
same systematic dependence on the number of fit param- 
eters. This is not surprising, as precise extractions of G^; 
are not limited to higher values. 

Finally, one can use the discrepancy between calcula- 
tions and measurements of the hyperfine splitting in the 
hydrogen ground state to extract the magnetic radius, as- 
suming that all of the other proton structure corrections, 
including the charge radius, are well known. The analysis 
by Volotka et al [95] yields a radius of 0.778(29) fm, in 
agreement with the Mainz result [40]. However, an up- 
dated calculation [45] which uses updated proton form 
factors and spin structure is also consistent with the 
hyperfine splitting even though it takes a larger mag- 
netic RMS radius of 0.858 fm, consistent with our re- 
sults. Thus, it appears that the hyperfine splitting does 
not have the sensitivity to precisely constrain the mag- 
netic radius, given the current uncertainty in the other 
portions of the proton structure corrections. 

So while the charge radius has a clear discrepancy be- 
tween muonic hydrogen lamb shift measurements and the 



various measurements using the electron-proton interac- 
tion, the situation for the magnetic radius is less clear. 
The extractions of the magnetic radius are more sensitive 
to small corrections, both for electron scattering and ex- 
tractions from hyperfine splitting measurements, making 
it difficult to tell at this point if there are, in fact, real 
discrepancies between the techniques. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we present an updated extraction of the 
form factor ratio i^pGe/Gm from the data of Ref. py]. 
We find a somewhat lower value for jupGe/Gm than the 
initial extraction for the entire dataset, consistent with 
two recent high precision measurements |4Q| , l4l| . The 
new analysis does not change our previous conclusion, 
i.e., that there is clear indication of a ratio smaller than 
unity, even for low Q^, indicating the necessity of in- 
cluding relativistic effects in any calculation of the form 
factors in this region. Both the form factors and pro- 
ton charge radius extraction from various e-p scattering 
measurements are in agreement, and there is still a sig- 
nificant disagreement with the charge radius as extracted 
from muonic hydrogen [43j . 
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